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The recent work by Achlioptas, D'Souza, and Spencer opened up the possibility of obtaining a 
discontinuous (explosive) percolation transition by changing the stochastic rule of bond occupation. 
Despite the active research on this subject, several questions still remain open about the leading 
mechanism and the properties of the system. We review the largest cluster and the Gaussian models 
recently introduced. We show that, to obtain a discontinuous transition it is solely necessary to 
control the size of the largest cluster, suppressing the growth of a cluster differing significantly, 
in size, from the average one. As expected for a discontinuous transition, a Gaussian cluster-size 
distribution and compact clusters are obtained. The surface of the clusters is fractal, with the same 
fractal dimension of the watershed line. 
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I. INTRODUCTION 

Percolation, the paradigm for random connectivity, has 
since Hammersley pQ been one of the most often applied 
statistical models [3]. Its phase transition being re- 
lated to magnetic models [3] is in all dimensions one of 
the most robust continuous transitions known. This ex- 
plains the enormous excitement generated by the recent 
work by Achlioptas, D'Souza, and Spencer [5] describ- 
ing a stochastic rule apparently yielding a discontinuous 
percolation transition on a fully connected graph. A dis- 
continuous percolation transition is observed when the 
growth of the largest cluster is systematically suppressed 
E] , promoting the formation of several large components 
that eventually merge in an explosive way [7]. Several 
aggregation models, based on percolation, have been de- 
veloped to achieve this change in the nature of the tran- 
sition El These models are generally classi- 
fied as explosive percolation, the name given in the orig- 
inal work that triggered the field [5]. In that work, a 
best-of-two product rule is proposed where, at each iter- 
ation, two unoccupied bonds are randomly selected but 
only the one which minimizes the product of the mass of 
the clusters, connected with the bond, is occupied. This 
work was originally studied by Achlioptas et al. [5] on 
a fully connected graph and analyzed in detail by Fried- 
man and Landsberg [7 . Ziff reported simulations on a 
regular square lattice [T3J [T3] , while Radicchi and For- 
tunato [13 [TB] and Cho et al. [T7] on scale-free net- 
works. However, reported results of finite-size studies 
and size distributions are not consistent with a discon- 
tinuous transition. For example, Ziff [T3], Radicchi and 
Fortunato [IB], as well as Cho et al. [5] found a power- 
law cluster-size distribution with an exponent close to 
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two. Although, different from the exponent of classical 
percolation, the sole fact of finding a power law is un- 
typical for discontinuous transitions. Also unusual for 
a discontinuous transition is that the clusters are frac- 
tal, as we found happens for the Achlioptas rule, from 
the behavior of the order parameter with the system size 
Q31 [TB] . Since then, various rules have been devised [BJ- 
[H HUj [TT] , all attempting a discontinuous transition to- 
wards an infinite cluster. In all proposed models one 
tries to keep the clusters of similar size and, for random 
graphs, the internal bonds of clusters should also be ad- 
ditionally suppressed El E] . For example, Moreira et al. 
[8] have proposed a Hamiltonian formalism which pro- 
vides a clear connection between equilibrium statistical 
mechanics and explosive percolation. They have shown 
that, for obtaining a discontinuous transition the size of 
the growing clusters should be kept approximately the 
same and, on random graphs, merging bonds (connect- 
ing different clusters) should dominate over the redun- 
dant ones (connecting sites in the same cluster). In this 
manuscript we review the work introduced by Araujo and 
Herrmann [Bj where two models are proposed yielding 
clear discontinuous transitions: the largest cluster and 
the Gaussian models. In the latter, at the percolation 
threshold, a bimodal cluster size distribution is found 
consistent with the nature of the transition. For both 
models, the cluster perimeters are fractal with a fractal 
dimension of 1.23 ± 0.03, similar to the one observed for 
watersheds PUDI] and other models [2"UH2"2] . 

More recently, the procedure proposed by Achlioptas 
et al. has been generalized to a best-of-m product rule in 
random graphs [531 H3] and regular lattices [35] to study 
its percolation and transport properties. The larger the 
set of bonds, m, the lower the probability that the occu- 
pied bond is related to the largest cluster, promoting the 
compactness of the percolation cluster, delaying the per- 
colation threshold, and above an intermediate value, im- 
proving the conductivity of the system [25] . Andrade et 
al. [2S] have shown that, at the percolation threshold, all 
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FIG. 1. Typical configuration of the system, at the percolation threshold, for the Largest cluster and the Gaussian models, 
both with a = 1. Pictures obtained from simulations on a square lattice with L 2 sites and L = 1024. Both models yield clear 
discontinuous transitions and their clusters are compact with the same fractal dimension of the watershed line. 



exponents for the size dependence of the spanning clus- 
ter, the conducting backbone, the cutting bonds, and the 
global conductance of the system, change continuously 
and significantly with m. 

A hybrid model has also been proposed [26| where an 
additional parameter is included to interpolate between 
the discontinuous transition, observed for m = 10, and 
the continuous one of classical percolation. The model, 
discloses a nonequilibrium tricritical percolation where 
explosive percolation is diluted with classical percolation. 
In the diagram for the model two transition lines were 
identified: a discontinuous and a critical line; both meet- 
ing at a tricritical point. In the work, the multicritical be- 
havior is characterized by a new set of critical exponents 
and a tricritical crossover between the discontinuous and 
the continuous regime is presented. 

Potential applications of explosive percolation are, for 
instance, the growth dynamics of the Human Protein Ho- 
mology Network |27j and the identification of communi- 
ties in real systems [2*5] . 

In this manuscript, we start with a description of the 
two models [largest cluster and Gaussian) in the next 
section. The nature of the transition and the fractal di- 
mension of the cluster perimeter is discussed in Sec. |III| 
with some final remarks in Sec. IIVI 



II. MODEL 

In the largest cluster model [6], bonds are randomly 
selected from the list of available ones. If, once occu- 
pied, the chosen bond would not lead to the formation 
or growth of the largest cluster, it is occupied; otherwise, 



the occupation occurs with probability, 



min < 1 , exp 



(1) 



where s is the size of the cluster obtained by occupying 
the bond, s is the average cluster size after the occupa- 
tion, and a is a parameter of the model that, for sim- 
plicity, we take equal to unity. With this parameter a it 
is possible to control the size distribution of the clusters. 
The larger the value of a, the lower the cluster-size dis- 
persion. For a < 0, all selected bonds are occupied and 
the model boils down to the classical bond percolation 
problem, with a continuous transition at the percolation 
threshold (see, for example, Ref. [2]). For a > 0, the for- 
mation of a largest cluster differing significantly in size 
from the average cluster size is systematically demoted, 
promoting the homogenization of the clusters size. This 
homogenization, induces the formation of a "powder keg" 
[55] which merges at the percolation threshold leading 
to a discontinuous transition. 

The Gaussian model |6j is generalization of the largest 
cluster one. While in the latter an occupation probabil- 
ity is solely defined to the bonds related with the largest 
cluster, and all the others are occupied with probability 
one, in the former, all bonds are occupied with a proba- 
bility given by Eq. ([I]) . At each iteration, an unoccupied 
bond is randomly selected and occupied with this prob- 
ability, which allows to explicitly control the cluster size 
distribution. We denote this model as Gaussian since 
the proposed expression corresponds to a Gaussian dis- 
tribution with average size s and size dispersion s/y/2. In 



3 




FIG. 2. Size dependence, at the percolation threshold, of the 
largest cluster model, of the susceptibility (x) and the order 
parameter (Poo) - fraction of sites belonging to the largest 
cluster. Results have been averaged over 10 4 realizations of 
square lattices with L 2 sites. The linear system size L ranges 
from 32 to 4096. The susceptibility scales linearly with the 
system size and the order parameter converges to a non-zero 
value in the thermodynamic limit. 

principle, any function constraining the formation of clus- 
ters differing significantly in size from the average clus- 
ter size could be considered leading to a discontinuous 
transition. Note that, in this model all clusters size are 
controlled while, in the largest cluster model solely the 
largest cluster is directly controlled and all the smaller 
ones can freely grow. 

Figure [T] has snapshots of typical configurations, at the 
percolation threshold, of both models. Obtained clusters 
are compact with fractal perimeters. In the next section, 
we discuss the discontinuous nature of the transition and 
the fractal dimension of the largest-cluster interface. 

III. RESULTS 

For nonequilibrium problems, where a free energy can- 
not be defined, transitions can still be classified based on 
the behavior of the order parameter 30 . A discontin- 
uous transition, is characterized by a jump in the order 
parameter, otherwise, a transition is denoted as contin- 
uous. For percolation, we define as order parameter the 
fraction of sites in the largest cluster (Poo) [2J. Here we 
also consider the second moment of the cluster size dis- 
tribution (x), defined as 

X=i£*?> ( 2 ) 

i 

where the sum runs over all clusters i. To estimate the 
percolation threshold we consider the average value of p 




FIG. 3. Size dependence, at the percolation threshold, of 
the Gaussian model, of the susceptibility (x) and the order 
parameter (Poo) - fraction of sites belonging to the largest 
cluster. Results have been averaged over 10 4 realizations of 
square lattices with L 2 sites. The linear system size L ranges 
from 32 to 4096. The susceptibility scales linearly with the 
system size and the order parameter converges to a non-zero 
value in the thermodynamic limit. 



(fraction of occupied bonds) at which a connected path 
linking opposite boundaries of the system is obtained. 
Considering different system sizes, for a = 1, we ob- 
tain for the percolation threshold, of the largest cluster 
model, p c = 0.632 ± 0.002. To identify the order of the 
transition, Fig. [2] presents a finite-size study for and 
X, averaged over 10 4 samples of square lattices with lin- 
ear sizes ranging from 32 to 4096. As we can see in the 
inset of Fig. [2j above a certain system size, the order 
parameter, at the percolation threshold, does not show 
any finite-size dependence, staying at a constant value in 
the thermodynamic limit (L — > oo). The second moment 
of the cluster size distribution (x) scales with L d (d = 2) 
which is a sign of a discontinuous transition [31] [32] • For 
the present models, the percolation thresholds are larger 
than the ones from previous models due to the compact- 
ness of the clusters (see Fig. [I]) . 

As example, for positive a, we present, in Fig. |3j a 
size dependence study of the order parameter and second 
moment of the cluster size distribution, for the Gaussian 
model, with a = 1, at the percolation threshold, on a 
regular square lattice with linear size (L) ranging from 
32 to 4096. Results were averaged over 10 4 samples. We 
extrapolate, for the infinite system, a percolation thresh- 
old p c = 0.56244 ± 0.00006. As for the largest cluster 
model, the density of the infinite cluster does not change 
significantly with the system size and the second moment 
of the cluster size distribution scales with L d (d = 2). As 
before, these results imply a discontinuous transition. 

As clearly seen in the snapshots of Fig. [T] clusters ob- 
tained with our models are compact but we find that the 
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FIG. 4. For the largest cluster and the Gaussian models 
(a — 1), the dependence on the stick size of the number of 
sticks necessary to follow the perimeter of the infinite cluster, 
with the yardstick method. Results have been averaged over 
10 4 realizations on a square lattice with 2048 2 sites. For vi- 
sual clarity, data for the Gaussian model have been vertically 
shifted by a factor of 0.1. 



FIG. 6. Cluster size distribution for the percolation thresh- 
old of the Gaussian model (a = 1), on a square lattice with 
L 2 sites. Results have been averaged over 10 4 samples and 
L = {1024,2048,4096}. To reduce finite-size effects, the con- 
tribution of the largest clusters has been neglected. Black- 
dashed lines are two Gaussian distributions fitting the re- 
sults from simulation. The black-solid line is the sum of both 
curves. 




FIG. 5. Snapshot of the watershed line obtained by randomly 
occupying all bonds in the system except the ones leading to 
an infinite connection, i.e., closing a path between the bottom 
and the top of the system. The fractal dimension of the inter- 
face line is the same as the one of the surface of the infinite 
cluster of the largest cluster and Gaussian models. 



surface is fractal. For the Gaussian model, we calculate 
for the cluster perimeter a fractal dimension of 1.23±0.03, 
obtained with the yardstick method [33] (Fig. 2|. For the 
largest cluster model, it is also characterized by a fractal 
perimeter with a fractal dimension of 1.26 ±0.04 (Fig. [4]). 
Compact clusters with fractal surface were also reported 
for irreversible aggregation growth in the limit of high 
concentration by Kolb et al. [31]. The value of this frac- 
tal dimension of percolation is intriguingly close to the 
one found for watersheds (1.211 ± 0.001) QHCES], ran- 
dom polymers in strongly disordered media (1.22 ±0.02) 
|20) . and several other models [35]. The simplest way to 
obtain the watershed line has been proposed by Cieplak 
et al. [35] and consists in randomly occupying bonds in 
the system by systematically suppressing the formation 
of a path connecting opposite borders of the system, i.e., 
any bond which leads to the formation of a cluster of 
connected sites touching the bottom and the top of the 
system is never occupied. As seen in Fig. [5] in the limit 
where all the other bonds are selected, only two clusters 
exist, separated by a watershed line. This line is fractal, 
with the same fractal dimension of the largest-cluster in- 
terface of the two models discussed here. 

Figure [6] shows the cluster size distribution, P(s,a), 
for different system sizes, obtained with the Gaussian 
model. Measurements have been performed at the per- 
colation threshold on a square lattice with 1024 2 , 2048 2 , 
and 4096 2 sites, and averaged over 10 4 samples. To re- 
duce finite-size effects, we neglect the contribution of the 
largest cluster. Two characteristic peaks are observed, a 
typical feature of a discontinuous transition. For a finite 
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system, at the percolation threshold of such transitions, 
coexistence of the percolative and non-percolative states 
is expected [57] . 

IV. FINAL REMARKS 

In this manuscript we summarize the main properties 
of two models of explosive percolation yielding clear dis- 
continuous transitions: the largest cluster and the Gaus- 
sian. With the largest cluster model we conclude that, 
on a regular lattice, to obtain an explosive transition is 
solely necessary to suppress the growth of the largest 
cluster what, indirectly, promotes the homogenization of 
the clusters size. With the Gaussian model the discon- 
tinuous nature of the transition is supported by a bi- 
modal cluster-size distribution resulting from the coex- 



istence of a percolative and non-percolative state. For 
both models, the fractal perimeter of the largest clus- 
ter is intriguingly close to the one found for watersheds 
(1.211 ± 0.001) [TS], random polymers in strongly dis- 
ordered media (1.22 ± 0.02) [50], optimal path cracking 
[211 [22] , and several other models [35] , and we have ar- 
guments that they are actually identical. 
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